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^J | Abstract. We develop a physical model of the zodiacal cloud incor- 

q^ ■ porating the real dust sources of asteroidal, cometary, and kuiperoidal 

OQ \ origin. Using the inferred distribution of the zodiacal dust, we compute 

its thermal emission and scattering at several wavelengths (1.25, 5, and 

20 /im) as a function of NGST location assumed to be at 1 AU or 3 AU. 

Areas on the sky with a minimum of zodiacal light are determined. 
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q< ' 1. Physical model for zodiacal thermal radiation and scattering 

on" 

Improvements in the zodiacal light emission and scattering to be given by space 
observations at 3 AU, compared with observations near the Earth, were dis- 
cussed by Mather and Beichman (1996). Unfortunately, a rather accurate multi- 
parametric model of the zodiacal brightness derived by Kelsall et al. (1998) from 
the COBE data cannot be reliably extrapolated to heliocentric distances as large 
as 3 AU. We have developed a physical model of the zodiacal cloud incorporating 
the real dust sources of asteroidal and cometary origin, which makes it possible 
to evaluate quantitatively the zodiacal light emission and scattering throughout 
/\ ' the Solar system (Gorkavyi et al. 1997a). This model considerably improves 

3 . our previous 'reference model' based on the use of the continuity equation for 

distribution function of dust particles (Gorkavyi et al. 1997b, c, 1998) and en- 
ables us to obtain more reliable results than a phenomenological modelling of 
the zodiacal light (e.g. Ebbets 1998). Below, we describe an improved model 
that represents a 3D-grid containing 45 x 180 x 244 = 2 • 10 6 cells with a step 
in (heliocentric latitude ip, longitude A, and radius R) to be (2°, 2°, 0.025 R 
[AU]). Using a new numerical approach to the dynamics of minor bodies and 
dust particles, we increase the number of particle positions employed in each 
model to 10 11 without using a supercomputer (Ozernoy et al. 2000, Gorkavyi 
et al. 2000). We compute here the distribution of the zodiacal dust emission 
and scattering in the Solar system. The processes influencing the dust particle 
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dynamics include gravitational scattering on all planets, except Pluto (inclina- 
tions and precession of the planets are neglected), mean motion resonances, and 
the Poynting-Robertson/solar wind drags. 

We employ here 931 sources of dust particles (with eccentricities e < 1), 
which include 284 asteroids, 451 short-period comets, and 196 Kuiper belt ob- 
jects. We adopt a two-component approximation for dust particle size distribu- 
tion, which is characterized by the parameter (3 = 0.285 and 0.057 for (l-2)/iin 
and (5-10)/xm particles, respectively. 

The inferred dust density distribution, n(R, <p, A), of the IPD at heliocentric 
distances 0.5 < R < 100 AU enables us to compute the scattered light and the 
thermal emission of the zodiacal cloud as a function of the observer's latitude 
and longitude. 

The thermal emission and scattering of the zodiacal cloud are given by: 

I v oc J n(R, <p, A) [(1 - A)B V {T) + AF®$(0)] ds 

where B v is the Planck function at frequency v, A is albedo, T is the dust 
temperature, T(R) = TqR- 5 with T = 286°K and 5 = 0.467 (Kelsall et al. 
1998), F® oc l/R 2 is the solar flux, <E>(0) is the phase function at scattering 
angle 6 (Kelsall et al. 1998). 

Note that although the phase function only weakly depends on v, the value 
of I u strongly depends on v through the solar flux. Besides, the Planck function 
strongly depends on heliocentric distance, especially in the near IR (1-5 jim). 

2. The zodiacal light at R=l and 3 AU 

Using the inferred distribution of the zodiacal dust, we have computed a variety 
of zodiacal light maps, both for thermal emission and scattered components, at 
different locations (R, Z) of the NGST assumed to be at 1 AU or 3 AU (see Fig. 1 
for the 5/um-emission of asteroidal dust and Fig. 2 for that of kuiperoidal dust). 
The brightness as a function of latitude ip and longitude A (in the telescope's 
frame) is given in logarithmic scale (the neighboring contour intensities differ 
by \J~e). The Sun's position is (0,0) for Z = and is shifted to a negative <p for 
Z = 0.25 AU. 

At each location, there is a minimum in the zodiacal light which can be seen 
as a 'dark spot' (or several spots) in the computed Figures. The positions of 
those minima are explained by an interplay between dependencies of emissivity 
upon density and temperature. 

Our numerical results concerning the computed zodiacal light at different 
locations of NGST are summarized in Tables 1 and 2. The data presented in 
Table 2 for R = 1 AU, 2nd column can be compared with the COBE data 
I [0°, 90°) /I (90°) = 3.5 (1.25/wn band), = 2.3 (5/zm band), and = 3.3 (20/xm 
band). Interestingly, it is the kuiperoidal component (hi) whose shape resembles 
most closely the profile of the zodiacal cloud, which might imply that the role of 
kuiperoidal particles has been highly underestimated so far. 



Table 1. Components of the model zodiacal cloud. 
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Table 2. Brightness of different components of the zodiacal cloud 

assuming the intensity of each component near the Earth be I((p = 90°, A) = 1 at the pole 
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"Type of dust listed in Table 1. 

"intensity at the north pole, I po i e 

'intensity at the ecliptic plane in the direction of 90° from the Sun, I(ip — 0°, A = 90°) 

intensity at the ecliptic plane in anti-sun direction, I{ip = 0°, A = 180°) 

d ^min(v 3 min, A = 180°), where <£>min is the latitude at which the zodiacal light brightness 
is minimal in the anti-sun direction (A = 180°) 

Vmin [°] 

■^All 5/im-intensities at R — 3 AU, Z — need to be multiplied by 10~ 3 



3. Conclusions 

1. The structure of the zodiacal dust cloud computed in the present work is 
substantially non-uniform: 
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(i) near the Earth, the thickness of the dust layer is the largest for cometary 
particles and the smallest for asteroidal particles, with kuiperoidal particles oc- 
cupying an intermediate position; 

(ii) the larger the size of asteroidal dust particles, the thinner is the layer 
comprised of such particles; for cometary particles a reverse, although a weaker, 
dependence takes place; 

(hi) the partial contribution of particles of different origin and size changes 
with heliocentric distance: as the distance increasess, the contribution of aster- 
oidal particles (especially of large size) into the zodiacal light emission sharply 
decreases, whereas the contribution of the kuiperoidal component increases; 

(iv) The latitudinal dependence of the zodiacal emission is different for 
different components of the dust cloud. In the anti-sun direction (A = 180°), 
the maximum of the asteroidal dust emission at 5 /im is reached at the ecliptic 
plane, whereas the brightness of the cometary component has its minimum here. 

2. Observations made from R = 1 AU, Z = 0.25 AU at all wavelengths 
would detect the brightness of the zodiacal light in the direction of the pole at 
level of 20-30% of that seen near the Earth. 

3. At R = 3 AU, the zodiacal light brightness in the direction of the pole 
decreases to 5 — 6% at 1.25 fim; to 0.05 — 0.07% at 5 /xm, and to 8 — 10% at 
20 /[/m compared to that seen near the Earth. 

4. At each wavelength, there are certain regions on the sky having a minimal 
zodiacal light brightness. At short wavelengths (~ 1.25 /um), these regions are 
around the poles, and at long wavelengths (such as 5 to 20 /im) these regions 
are shifted toward the ecliptic. This result needs to be taken into consideration 
while planning observations with NGST. 

5. Further improvements in our physical modelling of the zodiacal light are 
feasible. Using COBE and IRAS data, the contribution of dust components of 
different origin into the zodiacal cloud could be determined. Those improve- 
ments would make far going implications possible for extragalactic astronomy 
and cosmology. 

Acknowledgments. This work has been supported by NASA grant NAG5- 
7065 to George Mason University. N.G. acknowledges the NRC-NAS associate- 
ship. 

References 

Ebbets, D. 1998, "Zodical Light Model for NGST Studies", (unpublished) 
Gorkavyi, N.N., Ozernoy, L.M., & Mather, J.C. 1997a, ApJ 474, 496 
Gorkavyi, N.N., Ozernoy, L.M., Mather, J.C. & Taidakova, T.A. 1997b, BAAS 
29, 1310 

- 1997c, ApJ 488, 268 

1998, Earth, Planets and Space, 50, 539 

- 2000, Planetary and Space Science (submitted) 
Kelsall, T. et al. 1998, ApJ 508, 44 

Mather, J.C. & Beichman, C.A. 1996, "Unveiling the Cosmic Infrared Back- 
ground", Ed. E. Dwek, p. 271 
Ozernoy, L.M., Gorkavyi, N.N. & Taidakova, T.A. 2000, Icarus (submitted) 



ASTEROIDAL DUST, R=1 AU, Z = AU 
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Figure 1. Emission of a small (r = 1 — 2/xm) asteroidal dust in the 
5/um band 



KUIPEROIDAL DUST, R=1 AU, Z = AU 
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Figure 2. Emission of a small (r = 1 — 2/xm) kuiperoidal dust in the 
5/um band 
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